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1 INTRODUCTION 

Integrable models in two space-time dimensions provide a fascinating arena for investigating non- 
perturbative phenomena in quantum field theory |], |2], The integrable theories that we shall 
focus on here can be divided into three broad classes, but they all have in common the properties 
of asymptotic freedom and dynamical mass generation which make them akin to realistic models 
of particle interactions in four dimensions. First, there are bosonic sigma models based on a sym- 
metric space G/H. Any model of this sort is classically integrable, but unfortunately anomalies 
can destroy integrability at the quantum level unless H is simple. Second, there are fermionic 
models of Gross-Neveu or Thirring type with four-fermion interactions. Finally, there are hybrids 
of the two previous classes in which fermions are added to bosonic symmetric space models in 
various ways. In some cases the addition of fermions can cancel the anomalies encountered in 
the purely bosonic theory, leading to interesting new examples of quantum integrable theories; 
in particular this is believed to occur for supersymmetric sigma models based on a symmetric 
space G/H. A review of these matters with detailed references can be found in ||. 

Exact S-matrices have been proposed for many of the two-dimensional models which are 
thought to be integrable at the quantum level |3J, [IB], [TJ|, |T3|, [H], [T7|, [TJ|. These S-matrices 



describe the scattering of some conjectured set of particles states, and they are postulated on 
the basis of the specific symmetries of the model in question together with the usual axioms 
of S-matrix theory and the powerful constraint of factorization. (There may also be additional 
information available about the theory, such as the existence of bound-states.) However, such 



S-matrices are always subject to CDD ambiguities |19|, [Uj which cannot be further constrained 
by these general considerations. What is needed is some completely non-perturbative way of 
fixing the CDD ambiguities by testing the proposed equivalence between an S-matrix on the 
one hand and renormalized Lagrangian perturbation theory on the other. In this paper we shall 
explain such a programme which can be applied to any of the types of theory mentioned above. 

The technique, which was pioneered in ||, ^ Q following earlier work in |2y, [21], ^3|, p4f| , 
can be applied to any integrable model possessing a group G of global symmetries. The particle 
states in such a theory fall into representations of G, and the symmetry is generated by conserved 
charges in the Lie algebra. The idea is to couple the theory to some particular conserved charge 
Q by modifying the Hamiltonian from H to H — hQ, where h is a coupling constant of mass 
dimension 1 (so that Q is dimensionless) . The corresponding change in the ground-state energy 
density 5S(h) = S(h) — S(0) can be computed straightforwardly in perturbation theory. But 
in the special case of an integrable model, it can also be computed by a very different non- 
perturbative method starting from the exact S-matrix and using the Thermodynamic Bethe 
Ansatz (TBA). The comparison of these calculations provides a powerful check that the S- 
matrix and the Lagrangian really do correspond to one another and it also allows one to extract 
an exact expression for the mass gap of the model. We now describe more quantitatively how 
this comparison works. 

The result of a perturbative calculation of the ground-state energy density — which is of course 
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a renormalization group-invariant quantity — is an expansion 



SS(h) = h^Ujgih/Ky- 1 

3=0 

in ascending powers of the running coupling g(h/A) where the ctj are dimensionless numbers. 
The running coupling can be found by integrating the usual beta-function equation 

^=m = -t^ Obtain ^ = A ln£ + £ln1n£ + (lnln£/ln£), 

where the absence of a constant term in the solution defines the mass scale A for the particular 
renormalization scheme used. In the models we are considering, 0\ > and we expect pertur- 
bation theory to be valid in the asymptotic regime where \i ^> A. Combining the expressions 
above we find that the ground-state energy density is given by: 

S£(h) „ , h 02 , , h , h /, h\ 

— ¥- = a o 0i In - + a '-f In In - + an + CM In In -/ In - . 1 
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Notice that if ctQ 7^ 0> there is a classical or tree-level contribution to the ground-state energy 
density and this quantity is consequently unbounded as h becomes large. Otherwise, ao = and 
the leading contribution to the ground-state energy is then the constant term a±. We shall have 
more to say shortly about the important differences between these situations. 

Turning now to the TBA calculation: in order to extract the ground-state energy from the S- 
matrix, one must confront a set of coupled integral equations of Wiener-Hopf type. The special 
circumstances of interest to us are those in which the temperature is zero, with the coupling 
to the charge Q acting like a chemical potential. In general the resulting equations cannot be 
solved exactly, but it is possible, at least under certain simplifying assumptions, to generate an 
expansion of the solution in h/m valid when h^m, where m is some physical mass parameter 
occurring in the S-matrix. The result is 

= /t ln— + Kilnln— + k 2 + £>( In In — / In—), (2) 

where the Kj are dimensionless numbers. The exact values of the parameters Kj depend sensitively 
on the analytic structure of the S-matrix and, in particular, one finds that the presence or absence 
of CDD factors alters dramatically the result for the ground-state energy obtained in this way. 

Equality of the expressions (1) and (2) gives a powerful check that the Lagrangian and 
S-matrix descriptions are consistent. The greater the accuracy to which we can calculate these 
expressions, the more conclusive this check will be. Even finding just the first few terms, however, 
can provide strong evidence in favour of the proposed S-matrix used in the TBA calculation, 
allowing us to argue that any alteration by CDD factors would destroy the delicate agreement 
with the perturbative result. In addition, the consistency of (1) and (2) clearly determines the 
mass gap m/A, at least if we have calculated the expressions involved to sufficient accuracy. 
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At this stage it is useful to look more closely at the two possible cases to which we drew 
attention earlier. If a® 7^ 0, we see that the TBA calculation must reproduce the first two 
coefficients of the beta-function through the conditions: 

aoA = , ao/VA = K i > ( 3 ) 

providing a highly non-trivial test of the S-matrix. In these circumstances we can also read off 
the value of the mass gap in the model: 

ln(m/A) = («a - ai)/«o = («2 - a 1 )/a p 1 . (4) 

We emphasize that these relations are deduced by comparing the terms written explicitly in (1) 
and (2) and they therefore rely on just a one-loop perturbative calculation of the ground-state 
energy. If, instead, ao = 0, we must have Kq — Ki — and k 2 = a.\ for consistency and we are 
then unable to find the mass-gap by taking the expressions to the order given explicitly in (1) 
and (2). To extract m/A in this situation one needs to extend both the TBA and perturbative 
calculations to higher orders; in fact it is not hard to see that this would involve at least a 
three-loop perturbative calculation. 

The strategy we have just outlined has been applied to several series of integrable models, 
with results we shall summarize in section 3. The expressions for the exact mass-gaps derived 
in this way are very useful, because they provide bench-marks for the reliability of other non- 
perturbative approaches, such as lattice simulations. To minimize the amount of work involved, 
it is clearly desirable to try to choose Q so as to produce a classical term in the perturbative 
expression for the ground-state energy (1) because, as we have explained, the mass-gap can then 
be found from the TBA analysis in conjunction with a one-loop perturbative calculation. It 
can be shown |L2j that such a choice of Q is possible for any bosonic or supersymmetric sigma 



model based on a symmetric space G/H, a result which unifies the treatment of various examples 
considered previously in j|, [5|. |8|. [|. [1]], |TTJ . For the purely fermionic Gross-Neveu models, however, 
it seems that generically the classical contribution vanishes, and so a three-loop calculation is 
necessary in order to find the mass gap || [7j . 

There is another crucial consideration to be borne in mind when choosing Q. In each of 
the cases considered so far, an important simplifying assumption was made regarding the TBA 
calculation, namely, that for very particular choices of the charge Q, only a small number of 
particles — in fact those with the largest charge/mass ratio — contribute to the new ground-state. 
The original one-particle states can be chosen to be eigenvectors of the new Hamiltonian H — hQ 
with eigenvalues rrii cosh#j — hqi, where qi is the charge of the particle labelled by i and 6i is its 
rapidity. If i = 1 labels the particle type with the largest charge/mass ratio, then for sufficiently 
small values of h such that hqi < mi it is not energetically favourable to find any particles in 
the new ground-state, and so 5£{h) = 0. As A passes the threshold value mi/qi, it suddenly 
becomes favourable to fill the original vacuum with particles. For large h, it would seem that any 
arbitrary particle could appear in the new ground-state and hence the TBA calculation would 
have to keep track of all particles. However, for very particular choices of the charge Q it seems 
that the particles with largest charge/mass ratio actually repel other particles and so only they 
appear in the ground state. This assumption greatly simplifies the solution of the TBA equations 
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and in most of the existing papers it is taken as a working hypothesis which is vindicated by the 
consistency of the final results. With more care, it can actually be proven from an analysis of 
the full TBA equations ]T2|. 

Having introduced the general idea behind the technique, we shall discuss in the next sec- 
tion how the TBA equations can be derived under the simplifying assumption explained above. 
We shall then summarize the known results, concluding with a more detailed example which 
illustrates the most important points. 



2 THE TBA EQUATIONS 

In the cases where a single particle-type contributes to the ground state, the TBA analysis is 
rather straightforward. The elastic scattering of two particles of the same species is described 
by an S- matrix element which is just a phase S(8), where 9 is the rapidity difference of the 
incoming particles. Since particle number and momenta are conserved on interaction, it makes 
sense to consider single particle states. In the dilute regime, where the particles are on average 
well separated, the wavefunction of iV particles is built up from these single particle states and 
has the form 

ty(xi, . . . ,xn) = ®( x q)((Q) expim(xi sinh^i + . . . + x^ sinh^) 

QeS N 

where the sum is taken over all permutations Q = {Qx, ■ ■ ■ , Qn} of {1, . . . , iV} with 

q, x \ = { 1 if x Qi < X Q2 < ■ ■ ■ < x Qn 
I otherwise 

and where C(Q) are numbers defined via the S-matrix so that if Q differs from Q' only by an 
exchange of two specific elements of the list Qi and Qj, say, then 

C(Q') = S {Bi - OMiQ) . 

This means that the ({Q) are determined up to an overall factor. 
We impose periodic boundary conditions on the wavef unctions: 

^(xt, ...,Xj,..., x N ) = $?{xi, xj+L, x N ) 

which imply that, for each fixed j, 

exp(imL sinh 9j) S(9j — k ) = 1 . 

Taking the logarithm of these equations we have 

mL sinh 9j — i In S(9j — Ok) = 2ixn,j (5) 
Mi 
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where the {rij} are a set of N integers which we assume determine uniquely the set of rapidities 
{Oj} (this is guaranteed in the dilute regime where the sinh^- dominates in (5)). It is crucial that 
in all the integrable models we are considering here the particles behave as fermions in rapidity 
space. The condition (5) then proves to be the key ingredient when we come to analyse the 
thermodynamics of such a system. 

The thermodynamic limit can be taken in the usual way by letting N — > oo and L — > oo 
with N/L fixed. It now becomes meaningful to talk about the density of one-particle states in 
rapidity space, g(9), defined by dn = Lg{9)d9, and also the density of occupied states in rapidity 
space, cr(0). It is also convenient at this stage to define from the S-matrix the quantities K{9) 
and R{9) given by 

K(0)=6(0)-R(0) = ±:±\nS(8). (6) 
By differentiating the key condition (5) we deduce 

e {6) = — cosh9 + K*a(9) (7) 
2ir 

(where / * g{9) = S^dff f(9')g(9 - 9') for any two functions / and g). When K(9) = 0, (7) 
reduces to the usual expression for the density-of-states in a free theory. But when there is an 
interaction between the particles, (7) implies that the total density and the density of occupied 
states are coupled in a complicated manner. Of course for overall consistency we must have 
c(9) < q(0), which is guaranteed in the dilute regime since a (9) is then small compared with 
g(9). 

For a macrostate specified by the density cr(9), the value of H — hQ per unit length is 

d9 



/°° at) 
— (mcosh9-h)a(9). 
-oo 2% 



(where we assume the charge is normalized to unity on our preferred particle states). To find 
the ground-state energy density we must minimize S(a; h) with respect to a (9) subject to the 
constraint (7). The solution to this variational problem involves filling all the available states of 
lowest rapidity first, so a(9) = g(9) up to some Fermi rapidity ±9p. It can be shown that this 
problem can be solved in terms of an energy density e{9) which satisfies the TBA equation at 
zero temperature: 

e + {9) + R*e-{9) =mcosh9-h (8) 



with the notation 



_ , > 



no) 




otherwise. 



The solution to (8) is concave and negative between the values ±#p, which are fixed by the 
condition e(±#p) — 0. The final result for the way in which the ground-state energy density 
behaves as a function of h is then given simply by 

77? /*^F 

S£(h) = — d9 cosh9e(9) . (9) 

Z7T J-6 F 



5 



The problem of finding the ground-state energy density has thus been reduced to solving the 
Wiener-Hopf integral equation (8). It is not possible to find the solution exactly for arbitrary h. 
But we require the solution only in the asymptotic regime with h ^> m and it is explained in 
and || how to develop such an asymptotic expansion for the solution based on the original 
approach of P3| . The nature of the expansion depends crucially on whether or not the Fourier 
transform R(u) of the kernel R(0) defined via 

r°° dm 

Rift) = / —cob(u6)R(u), 

J0 7T 

vanishes at the origin. Suppose that -R(O) = 0, and that we can decompose 1/ R{u) = G + {uS)G-{uS) 
where G±(u) are analytic in the upper/lower half planes with G+(u) = G-(-uj). It can be shown 
that if G+(i£) has an expansion for small £ like 

G+ft) = A exp(-a£lnO (l - &£ + 0(£ 2 )) (10) 

then the ground-state energy density for h 3> m takes the form (2) with 

Ko = Kt/Ko = 0+1/2 , k 2 /k = ln(v / 2^/G+(i)) - l + lnfc + a(7 E -l + ln8)-6 . (11) 

This matches the expansion in perturbation theory with a classical term in (1). If R(0) ^ 0, on 
the other hand, then the asymptotic expansion proceeds slightly differently and the result is an 
expression of the form (2) with k,q = K\ = 0. This matches the perturbative result (1) in the 
case when there is no tree-level contribution. 

At this point we should remark that in the regime h ^> m considered above, the Fermi 
rapidity is large and so the ground-state will contain a large number of particles. It would seem 
therefore that the system is far from the dilute situation used in our naive derivation of the TBA 
equations. It is one of the miracles of the TBA that it seems, nevertheless, to be valid for many 
systems in the deep ultra-violet regime, although the exact reason for this is not understood. 

The simple TBA analysis which we have presented also assumed that there is only one particle 
type that contributes to the new ground-state. In general many particles with different quantum 
numbers and masses can contribute and the TBA analysis then becomes much more complicated. 
If the scattering of the particles is purely elastic, the one-particle analysis can be generalized in 
an obvious way: there is a function e a (9) for each particle type with mass m a , and the TBA 
system (8) becomes a matrix equation involving 

where S a t,(8) is the S-matrix element between particles a and b. If the scattering is not elastic, 
there are still greater complications, with the net result being that the TBA system involves 
additional "magnon" degrees of freedom which behave like particles with zero mass. Nevertheless 
the analysis can sometimes be carried out successfully in this situation too, as we shall see later 
for a particular example. 
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3 SUMMARY OF RESULTS 



We now summarize the results obtained for several families of integrable models. In each case we 
shall define the theory by a Lagrangian £ displaying a global symmetry group G or GxG, where 
G = SU(n), SO(n) or Sp(n). To express the results compactly, we introduce the quantity 1/A 
= n, n—2 or 2n+2 respectively (this is the dual Coxeter number for SU(n) and SO(n) but twice 
the dual Coxeter number for Sp(n)). The fields in the Lagrangian transform in the defining 
representation of the symmetry group unless we state otherwise. The mass-gap m/A will be 
given for particles which also belong to the defining representation of the symmetry group. The 
renormalization scheme is MS except in one case. 

(i) O(n) sigma model || |j: 

2^ 

where (f) a is an n-component real scalar field obeying a a = 1- The TBA calculation confirms 
the S-matrix proposed in [Fj|, correctly predicting fa = l/27rA and fa = l/47r 2 A. The mass gap 



(8/e) A 

m ~ r(i + A) MS 



is consistent with the 1/n expansion |26] 
(ii) GxG principal chiral model §| | 



£ = - Tr (dJJd^U- 1 
9 v 



where U is a G-valued field. The TBA calculation with the S-matrices proposed in |T6| correctly 
predicts fa = l/16irxA, fa = 01/2 where x is the Dynkin index of the defining representation of 
G (x = 1/2 for SU(n), Sp(n); x = 1 for SO(n)). The mass gap is 

m = 2 (,A + i/ 2) ^A) A 



where d is the dimension of the defining representation of G (d — n for SU(n), SO(n); d = 2n 
for Sp(n)). The mass gap has been measured on the lattice for G = SU(3) and the agreement 



with the theoretical value is quite accurate |p0 |. 
(iii) O(n) Gross- Neveu model (n > 4) H: 



where ip a is an n-component Majorana spinor. The S-matrices proposed in [y| [TJ| are found to 
be consistent with the TBA/perturbation calculation which correctly predicts fa = l/2nA and 
fa = — l/47r 2 A. The mass gap is 

(2e) A 
m ~ r(l-A) MS 
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which is consistent with the 1/n expansion |[27| . Note that our coupling constant g and mass-scale 
A^jg differ slightly from those used in |6|, [27| . 



(iv) SU(n) chiral Gross- Neveu model Jj]]: 
where ip a is an n-component complex spinor. The S-matrix for the massive sector conjectured 



in 



and (3\ 



15 ] is consistent with the TB A/perturbation calculation which correctly predicts f3\ = 1/n A 
1/n 2 A. The mass gap 

(e/4) A / 2 



m 



A 



r(i-A) 

agrees with the 1/n expansion J7|. See for full details of the renormalization scheme and note 
also our slightly different definitions of g and A. 

(v) O(n) super symmetric sigma model (n > 4) |T0[ : 

1 



C 



1 

29 



where the fields <p a and ip a are an n-component real scalar and spinor satisfying the constraints 
4>a4>a = 1 and (fiaipa = 0. The model is invariant under N = 1 supersymmetry transforma- 
tions mixing the bosons and fermions. The S-matrix conjectured in |17] is consistent with the 
TB A/perturbation calculation which correctly predicts (3\ = l/27rA and 02 = 0. The mass gap 



m 



2 2A sin(TrA) A 

VTA 



MS 



is consistent with the 1/n expansion |2 

(vi) SU(n) or CP™ -1 supersymmetric sigma model JIT ]: 

1 



1 
2~9 



4 



where = ^(z*d^z a — Zad^z*) and the fields z a and ip a are an n-component complex scalar and 
Dirac spinor satisfying the constraints z*z a = 1 and z*ip a = 0. The Lagrangian has a local U(l) 
invariance under which all complex fields transform by a phase and A^ transforms as a gauge 



field. The model is also invariant under N = 2 supersymmetry. The S-matrix proposed in JT8 
correctly predicts /3 X = 1/n A and /3 2 = 0. The formula for the mass-gap is 



sin(TrA) 
ttA ms 



which is consistent with the 1/n expansion |29| and with other non-perturbative approaches 



specific to N = 2 supersymmetric theories [31 
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4 DETAILED EXAMPLE: THE SUPERSYMMETRIC 
CP 71 ' 1 MODEL 

Most of the features of the general method we have outlined here are nicely illustrated by the 
last example, the supersymmetric CP n_1 model [32], so we now discuss this case in more 
detail. 

To introduce the S-matrix, we assume that in the quantum theory there exist states \a,i,6) 
representing fundamental particles and states \d,i,9) representing fundamental antiparticles. 
Here 9 is rapidity; % = 0, 1 distinguishs "bosons" and "fermions" (actually these particles carry 
fractional statistics |18|, [361); and indices a and a label the n and n representations of SU(n). 
The fundamental anti-particles can be regarded as bound-states of the fundamental particles, or 
vice- versa — an example of "nuclear democracy" . There are additional bound states transforming 
in all the antisymmetric representations of SU(n). 



The integrability of the model [3~3] implies that the S-matrix factorizes and that all S-matrix 
elements can be deduced from the two-body ones; furthermore, the S-matrix for any desired set 
of particles can be obtained from the S-matrix for the fundamental particles. Thus the entire 
S-matrix is specified by the amplitude 

(c, k, 9 2 ; d, 1, 9 U out|a, i, 9 X ; 6, j, 9 2 , in) = ^ =2 (^ - 9 2 )^S CGN (9 1 - 9 2 ) c a d b (12) 



which was first proposed by Koberle and Kurak [18]. Following [[35, |36|, we have written this 
proposal in a factorized form in which Scgn is the S-matrix for the chiral Gross-Neveu model 
T5| , which specifies the scattering of the SU(n) degrees of freedom, and Sn =2 controls the N = 2 



supersymmetric degrees of freedom, as described in [[36[] (although we have chosen to include the 
physical strip pole in the CGN part of the S-matrix). Explicit expressions for these factors are 



given in [11] 



The conjectured S-matrix is minimal in the sense that it has the minimum number of poles 
and zeros on the physical strip (the region < lm(9) < 7r) consistent with the requirements 
of symmetry, the existence of a bound-state and the axioms of S-matrix theory. But this still 
leaves open the possibility of adding CDD factors to the S-matrix [TP], |HJ; these spoil none 
of the axioms, they introduce no new poles on the physical strip and they passively respect the 
bootstrap equations. For our model the CDD ambiguities correspond to multiplying the S-matrix 
of the fundamental particles by factors of the form 



si** (§-&«) sinh (f-fK 2 
sinh (f + f«) sinh (f + K 2 



a. 



a: 



(13) 



where < a < 2. One of the conclusions of our analysis will be that the minimal form is the 
true S-matrix of the theory, so that all CDD factors are ruled out. 

We have emphasized that in following the general method for testing the equivalence between 
an S-matrix and a Lagrangian, there are at least two important points to be borne in mind when 
choosing the coupling to a conserved charge Q. First, in order to simplify the analysis of the 
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TBA system we must choose Q so that the new ground-state consists of a restricted number of 
particle types. In the models considered in [§, || || (7), §, |{| it was possible to find a generator Q 
such that the new ground-state contained a single particle type. In the present theory, however, 
we know that the lowest energy states must come in degenerate supersymmetric multiplets, since 
supersymmetry commutes with the SU(n) invariance. We might be tempted to make the same 
choice for Q as in the SU(n) principle chiral model || and the chiral Gross- Neveu model 0: 

Q = diag(l, (14) 

V n — 1 n—1/ 

for which there will be a single fundamental doublet \l,j,9) with the largest charge/mass ratio. 
But it turns out that this violates our second criterion regarding the choice of Q, which is that 
there should be a classical or tree-level term in the ground-state energy so that we can extract 
the mass-gap via a one-loop, rather than a higher-loop, calculation. We are therefore motivated 
to consider an alternative choice: 

Q = diag(l,-1,0,...,0). (15) 

This does indeed lead to the desired tree-level term, but the TBA analysis is complicated by the 
fact that there are now two fundamental doublets with the largest charge/mass ratio, namely 
\l,j,9) and \2,j, 9), Despite this complication, the calculation proves tractable with this choice, 
as we shall describe below. 



4.1 Perturbation theory calculation 

The coupling of the theory to the charge (15) by a change in the Hamiltonian H — > H — hQ can 
be achieved by making a corresponding replacement do — > <9o + ihQ in the Lagrangian. Since 
we are interested in performing a one-loop calculation of the change in the ground-state energy 
density as a function of h, it is enough to expand the resulting Lagrangian to quadratic order 
in an independent set of fields, and we can drop all terms which are independent of h to this 
order. By exploiting the local U(l) invariance of the Lagrangian, we can take z\ to be real 
and we can solve the bosonic constraint z a z* = 1 by writing Z\ = {(1 — |vr| 2 )(| + 0)} 1 / 2 and 
Z2 = e %e {{\ — |vr| 2 )(| — 0)} 1 / 2 where ir = (23, . . . , z n ) and 9, (j) are real. The fermionic degrees of 
freedom and the variable 9 decouple to quadratic order and we are left with the expression 

A-loo P = ^ {{d,<P? + l^vrl 2 + h 2 - 4^ 2 2 - h 2 \n\ 2 } . 

We see that there is indeed a tree-level term, as desired. 

Using standard dimensional regularization with the MS-scheme gives a one-loop expression 
for the ground-state energy 

h 2 h 2 h 2 n r /i 2 l 

SS{h) = ~ --ln2 + ^ 1-ln^ . 

2g 7T 47T [ jjL Z 
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It is important that when we substitute for the running coupling with the known values of the 
beta-function coefficients 

Pi = n/n, P 2 = , (16) 

the /i-dependence must cancel (to leading order) since the quantity we are computing is renormal- 
ization group-invariant. From above, the values of the other dimensionless numbers appearing 
in (1) are, for this model, 



a = -1/2 , oti = [n/iir) - (1/tt) In 2 . 
The explicit expression for the ground-state energy is thus, to the required order, 



;i7) 



58(h) 



h 2 n 
'~2T 



In 



h 



1 2 

- + - In 2 

2 n 



4.2 TBA calculation 



We must now consider the TBA equations for the model and solve them in the limit h ^> m. 
Following the hypothesis introduced earlier, we assume that only the multiplets 9) and 
|2,j, 9) contribute to the ground-state. Since the scattering of these multiplets is purely elastic, 
it is not necessary to perform a diagonalization in the space of SU(n) quantum numbers (although 
this diagonalization can be done |]25|| ). The remaining difficulty is that the S-matrix for these 
favoured states is still non-diagonal in the supersymmetric subspace. Fortunately, this problem 
can be solved — in fact it has been shown by Fendley and Intriligator ]36] that it is equivalent to 
diagonalizing the transfer matrix of the six vertex model at the free fermion point. 

The TBA equations involve energy densities e\(6) and e^(9) for the two supermultiplets and 
two magnon energy densities £o(#) and £q(9) which reflect the non-elastic nature of the scattering 
amongst the supersymmetric degrees of freedom (two magnons because we are dealing with N = 2 
super symmetry). The resulting TBA equations are, at T = 0, 



e a {0) - <pab * e b (9) - (p a i * £j (9) = m cosh 9 - h, 
m-<f>ai*e-(9) =0, 

where a = 1,2 and I = 0, 0. In terms of these the ground-state energy density is 



58(h) 

The kernels appearing in (19) are 

<M#) 



2tt 



(19 



eT(9) + e^(9) 



cosh^. 



(19) 



(20) 



1 d 



2m d9 



InScGN^ab) 



1 



sin(7r/ri) 



'20 



(0) 



>10 



(0) 



2n cosh 9 — cos(7r / n) 
1 sin(7r/n) 
27r cosh 9 + cos (tt / n) 
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where ScGN(@) b a i are the appropriate S-matrix elements of the chiral Gross-Neveu model. 

To simplify (19) it is important to notice that <t> a i{9) is a positive kernel, which implies that 
the magnon variables are given by ^(9) = and £f(t) = 4> a i * e~ (9). Furthermore, the solution 
does not distinguish between the values of the favoured SU(n) quantum numbers and so we have 
ei(9) = €2(9) = e(9). The four equations in (19) then reduce to a single equation for e{9) of the 
form (8) with kernel 

R{9) = 5(9) - (j> u (e) - MB) - [<f> w + 1O ] * [0io + 0io] (9) , 

and, from (20), each of the densities ei and eg contributes an amount to the ground-state energy 
given by (9). We have therefore succeeded in reducing the problem to the relatively simple case 
which we already know how to handle. 

Following the subsequent steps outlined in section 2, we find the Fourier transform of the 
kernel: 

- . , cosh((i — sinh(-7ra;) , 

H{u) - -j- exp ^-ku) . 

cosh (^nuj) 

This vanishes at the origin and we can decompose it as 1 / '(G + (u)G -{uj)) where G±(u) are 
analytic in the upper/lower half planes and G^(uj) = G + (—uj). The unique solution is 

G+(w) = T{ * - %~^ { \ - exp { I Inn - 

+ i.[h2 + i + (H)ln(H) + >i]} 
and we can indeed find an expansion of G+(i£) for small £ of the form (10) with 

k = ^ ) a = -1/2, 6 = l(l- 7E )-Iln4n+(|-i)ln(|-I) . 
Substituting in (11) and adding the contributions for the densities e\ and eg produces 

K = -n/2ir, «i = 0, k 2 /ko = ln(ra/7r) + lnsin(7r/n) + (21n2)/n- 1/2, (21) 
or explicitly, to the required order, 



, h 1 2 1 o 1 f n \ 1 • f 71 

In h - In 2 + In - + In sin - 

m 2 n \irj \n 



(22) 



4.3 Comparison of calculations 



By substituting the values (16), (17) and (21) into the general relations (3) and (4) — or simply 
by comparing the explicit final results (18) and (22) — we see that the TBA calculation correctly 
reproduces the universal coefficients of the beta-function and that it predicts the value of the 
mass-gap for the supersymmetric CP n_1 model to be m/A^g = (n/ir) sin(7r/n), as claimed. 
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It is instructive to consider what would have happened if we had carried out our analysis 
using the other charge (14), rather than (15). In that case only the multiplet \l,j,9) appears in 
the ground-state and the resulting TBA equations are simpler in as much as they involve only 
this single doublet, rather than two doublets. The system can be reduced to a single integral 
equation in a similar way, but on doing so we find a different kernel: 

R(9) = 5(9) - n (0) - 1O * 4> w {9) - lO * l5 (0) ■ 

The Fourier transform of this kernel does not vanish at the origin, and so, as explained previously, 
we would need to go beyond one-loop perturbation theory to carry out a non-trivial test of 
the S-matrix. This matches precisely the fact that with this different choice of charge the 
perturbative expansion of the ground-state energy density is also markedly different with no 
tree-level contribution. 

It is also interesting to see for this particular example how the calculation resolves the problem 
of CDD ambiguities in the S-matrix. An additional CDD factor of the form (13) would alter the 
kernel appearing in the TBA equation from R(uj) to 

6/ x o cosh((|-^)7r^) cosh(^(a - 1)W) 

K[U) - 2 —| r . 

cosh (^tt^J 

But this expression fails to vanish at the origin and so the argreement with the perturbative 
result is destroyed. 
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